Shingle 2.0: generalising self-consistent and automated domain discretisation for multi-scale geophysical models

The approaches taken to describe and develop spatial discretisations of the domains required for geophysical simulation models are commonly ad hoc, model or application specific and under-documented. This is particularly acute for simulation models that are flexible in their use of multi-scale, anisotropic, fully unstructured meshes where a relatively large number of heterogeneous parameters are required to constrain their full description. As a consequence, it can be difficult to reproduce simulations, ensure a provenance in model data handling and initialisation, and a challenge to conduct model intercomparisons rigorously. This paper takes a novel approach to spatial discretisation, considering it much like a numerical simulation model problem of its own. It introduces a generalised, extensible, self-documenting approach to carefully describe, and necessarily fully, the constraints over the heterogeneous parameter space that determine how a domain is spatially discretised. This additionally provides a method to accurately record these constraints, using high-level natural language based abstractions, that enables full accounts of provenance, sharing and distribution. Together with this description, a generalised consistent approach to unstructured mesh generation for geophysical models is developed, that is automated, robust and repeatable, quick-to-draft, rigorously verified and consistent to the source data throughout. This interprets the description above to execute a self-consistent spatial discretisation process, which is automatically validated to expected discrete characteristics and metrics.Comment: 18 pages, 10 figures, 1 table. Submitted for publication and under revie

Ocean model resolution dependence of Caribbean sea-level projections

Abstract Sea-level rise poses severe threats to coastal and low-lying regions around the world, by exacerbating coastal erosion and flooding. Adequate sea-level projections over the next decades are important for both decision making and for the development of successful adaptation strategies in these coastal and low-lying regions to climate change. Ocean components of climate models used in the most recent sea-level projections do not explicitly resolve ocean mesoscale processes. Only a few effects of these mesoscale processes are represented in these models, which leads to errors in the simulated properties of the ocean circulation that affect sea-level projections. Using the Caribbean Sea as an example region, we demonstrate a strong dependence of future sea-level change on ocean model resolution in simulations with a global climate model. The results indicate that, at least for the Caribbean Sea, adequate regional projections of sea-level change can only be obtained with ocean models which capture mesoscale processes.info:eu-repo/semantics/publishe

Water motion and vegetation control the pH dynamics in seagrass-dominated bays

info:eu-repo/semantics/publishe

On accurate momentum advection scheme for a z-level coordinate models

Abstract In this paper, we focus on a conservative momentum advection discretisation in the presence of zlayers. While in the 2D case conservation of momentum is achieved automatically for an Eulerian advection scheme, special attention is required in the multi-layer case. We show here that an artificial vertical structure of the flow can be introduced solely by the presence of the z-layers, which we refer to as the staircase problem. To avoid this staircase problem, the z-layers have to be remapped in a specific way. The remapping procedure also deals with the case of an uneven number of layers adjacent to a column side, thus allowing one to simulate flooding and drying phenomena in a 3D model

Comparison of unstructured, staggered grid methods for the shallow water equations

Unstructured grid models are receiving increased attention mainly because of their ability to provide a flexible spatial discretization. Hence, some areas can be resolved in great detail while not over-resolving other areas. Development of these models is an ongoing process with significant longstanding issues with spurious computational modes, efficiency, advection and Coriolis approximations, and so forth. However, many of these problems have been solved with the current generation of models which have much promise for coastal to global scale ocean modelling. Our purpose is to intercompare a class of unstructured grid models where the continuity equation reduces to a finite volume approximation. The momentum equations can be approximated with finite difference, finite element, or finite volume methods. Each of these methods can have advantages and disadvantages in different classes of problems that range from hydraulics to coastal and global ocean flows. Some of the more important differences are restrictions on grid irregularity and stability of the Coriolis term. The finite element version of the model has important advantages in the discretization of the Coriolis term and does not require a reconstruction of a tangential velocity component. The comparison is illustrated with a simple test case. (C) 2008 Elsevier Ltd. All rights reserved

Spatial Variations of Antarctic Intermediate Water in the Caribbean Sea Due To Vertical Mixing Along Its Path

Because of its pronounced fresh signature, the properties of the northward-flowing Antarctic Intermediate Water (AAIW) affect the Atlantic Meridional Overturning Circulation. Hence, understanding modifications of AAIW along its path is important. Here, we analyze AAIW changes along its path in the Caribbean Sea and assess whether vertical fluxes from background turbulence and from double-diffusive mixing in thermohaline staircases can explain these variations. We deduce the occurrence rate of staircases (7%) and estimate the flux ratio (Formula presented.) from Argo float profiles. In combination with vertical fluxes from background turbulence, these values are used in a steady-state advection-diffusion model to estimate the effective diffusivity of salt that arises from double diffusion (Formula presented.). This value for (Formula presented.) is similar to observed values (Schmitt, 2005, https://doi.org/10.1126/science.1108678), implying the observed modification of AAIW in the Caribbean Sea may be attributable primarily to vertical mixing in the region itself